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Electronic band structure and optical properties of zinc monochalcogenides with zinc-blende- and 
wurtzite-type structures were studied using the ab initio density functional method within the LDA, 
GGA, and LDA-t-C/ approaches. Calculations of the optical spectra have been performed for the 
energy range 0-20 eV, with and without including spin-orbit coupling. Reflectivity, absorption 
and extinction coefficients, and refractive index have been computed from the imaginary part of 
the dielectric function using the Kramers-Kronig transformations. A rigid shift of the calculated 
optical spectra is found to provide a good first approximation to reproduce experimental observations 
for almost all the zinc monochalcogenide phases considered. By inspection of the calculated and 
experimentally determined band-gap values for the zinc monochalcogenide series, the band gap of 
ZnO with zinc-blende structure has been estimated. 

PACS numbers: 71.15.-m; 71.22.-|-i 

Keywords: Zinc monochalcogenides, electronic structure and optical spectra 



I. INTRODUCTION 

The zinc monochalcogenides (ZnA; A = O, S, Se, and 
Te) are the prototype II- VI semiconductors. These com- 
pounds are reported to crystallize in the zinc-blende(-z) 
and wurtzite (w) type structures. The ZnA-z phases 
are optically isotropic, while the ZnX-w phases are 
anisotropic with c as the polar axis. ZnA phases are 
primary candidate for optical device technology such as 
visual displays, high-density optical memories, transpar- 
ent conductors, solid-state laser devices, photodetectors, 
solar cells etc. So, knowledge about optical properties of 
these materials is especially important in the design and 
analysis of ZnA-based optoelectronic devices. 

Optical parameters for some of the ZnA phases have 
widely been studied experimentally in the past. De- 
tailed information on this subject is available for ZnO- 
^^^^^^^i^^ ZnS-w^ ZnS-z^aodi ZuSe-z^iS and 
ZnTe-Z)^ d°d^d3 ggg ^j^g systematized survey in Ref. [T^ . 
However, there are no experimental data on optical prop- 
erties of ZnSe-w, ZnTe-w, and ZnO-z. Furthermore, 
there is lack of consistency between some of the exper- 
imental values for the optical spectra. This is demon- 
strated in Fig. [U which displays reflectivity spectra for 
ZnO-w measured at T=300 K by three different groups. 
Dielectric response functions were calculated using the 
Kramers-Kronig relation. As is seen in Fig. [1] inten- 
sity of the imaginary part of the dielectric function (62) 
and reflectivity {R) corresponding to the fundamental 
absorption edge of ZnO-w are higher— than those at the 
energy range 10-15 eV, while in Ref. [ij it is vice verse. 
The optical spectra in Fig. [T] measured using the linearly 
polarized incident light for electric field {E) parallel (||) 
and perpendicular (±) to the c-axes are somehow close 
to those of Ref. [t! using non-polarized incident light. 

Using the experimental reflectivity data, a full set of 



optical spectra for ZnO have been calculatecP^ for the 
wide energy range 0-26 eV. Density-functional theory 
(DFT)i^ in the local-density approximation (LDA)^^ has 
also been used to calculate optical spectra for ZnO-w— 
and ZnS- by linear combination of atomic orbitals, 
and for ZnS-a^^ and ZnSe-ai^ by self-consistent linear 
combination of Gaussian orbitals. The optical spectra of 
ZnO (including excitons) has been investigated^ by solv- 
ing the Bethe-Salpeter equation. Band-structure studies 
have been performed by linearized augmented plane- wave 
method plus local orbitals (LAPW-I-LO) within the gen- 
eralized gradient and LDA with the multiorbital mean- 
field Hubbard potential (LDA-f [/) approximations. The 
latter approximation is found to correct not only the en- 
ergy location of the Zn-Sd electrons and associated band 
parameters (see also Refs . [21I and [22|) . but also to improve 
the optical response. Despite the shortcoming of DFT in 
relation to underestimation of band gaps, the locations 
of the major peaks in the calculated energy dependence 
of the optical spectra are found to be in good agreement 
with experimental data. 

It should be noted that the error in calculation of the 
band gap by DFT within LDA and generalized-gradient 
approximation (GGA) is more severe in semiconductors 
with strong Coulomb correlation effects than in other 
solids This is due to the mean-field character 

of the Kohn-Sham equations and the poor description 
of the strong Coulomb correlation and exchange inter- 
action between electrons in narrow d bands (viz. the 
potential U). Not only the band gap {Eg), but also the 
crystal-field (CF) and spin-orbit (SO) splitting energies 
(AcF and Aso)i the order of states at the top of the 
valence band (VB), the location of the Zn-3c? band and 
its width, and the band dispersion are foun d^^i^^'^^i^^ 
to be incorrect for ZnO-w by the ab initio full poten- 
tial (FP) and atomic-sphere-approximation (ASA) lin- 
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FIG. 1: Reflectivity spectra R{lo) for ZnO-w determined experimentally at 300 K in Refs. I9lll3 (solid circles), Ref. (open 
circles), and Ref. ,7 (solid lines) along with imaginary part of the dielectric response function (e2(<^)) calculated using the 
Kramers-Kronig relation. The results of Ref. [Tj (open circles) are used for both -E || c and iJ _L c, because no polarized incident 
light was used in the experiments. 



ear muffin-tin orbital (LMTO) metliods within the pure 
LDA,^^'^^ and by the projector-augmented wave (PAW) 
method within LDA and GGAj ^^i^^ These findings were 
ascribed^'^'^^ to strong Coulomb correlation effects. DFT 
calculations within LDA plus self-interaction correction 
(LDA+5/C) and LDA+U are iouiiSL2^ to rectify the 
errors related to Acf and Aso, order of states at the top 
VB, width and location of the Zn 3d band, as well as 
effective masses. In other semiconductors, in which the 
Coulomb correlation is not sufficiently strong, the Acf 
and Ago values derived from DFT calculations within 
LDA are found to be quite accurate. This was demon- 
strated for diamond-like group IV, z-type group III-V, 
II-VI, and I- VII semiconductors w-type AIN, GaN, 
and InN2£ using the LAPW and VASP-PAW, the w- 
type CdS and CdSe^ z-type ZnSe, CdTe, HgTe^ us- 
ing the ab initio LMTO-ASA, z- and w-type ZnSe and 
ZnTe?l''?? as well as z-type CdTe^^ using the VASP-PAW 
and FP LMTO methods. Ahhough the SO sphtting at 
the top of VB is known to play an important role in 
electronic structure, and chemical bonding of semicon- 
ductor s2ii2^iS&i2^i2Si^iHiS there is no systematic study 
of the role of the SO coupling in optical properties of 
these materials. 

Several attempts have been undertaken to resolve the 



DFT eigenvalue problem. One such approach is uti- 
lization of the GW approximation ("G" stands for one- 
particle Green's function as derived from many-body per- 
turbation theory and "W" for Coulomb screened inter- 
actions) . Although GW removes most of the problems of 
LDA with regard to excited state properties, it fails to 
describe the semiconductors with strong Coulomb corre- 
lation effects. For example, two studies of the band gap of 
ZnO calculated using the GW correction underestimated 
Eg by 1.2 eV^ and overestimated it by 0.84 eVi^ Calcu- 
lations for Zn, Cd, and Hg monochalcogenides by the GW 
approach showed^ that the band-gap underestimation is 
in the range 0.3-0.6 eV. Combination of exact-exchange 
(EXX) DFT calculations and the optimized-effective GW 
potential approach is found"^^ to improve the agreement 
with the experimental band gaps and Zn-3d energy lev- 
els. Band gaps calculated within the EXX treatment are 
found to be in good agreement with experiment for the 
s-p semiconductors. Excellent agreement with exper- 
imental data was obtained^ also for locations of energy 
levels of the d bands of a number of semiconductors and 
insulators such as Ge, GaAs, CdS, Si, ZnS, C, BN, Ne, 
Ar, Kr, and Xe. 

Another means to correct the DFT eigenvalue error is 
to use the screened-exchange LDA.'"^ Compared to LDA 
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and GW, this approximation is found to be computa- 
tionally much less demanding, permitting self-consistent 
determination of the ground-state properties, and giving 
more correct band gaps and optical properties. Other 
considered approaches for ab initio computations of op- 
tical properties involve electron-hole interaction,"^^ par- 
tial inclusion of dynamical vertex corrections that neglect 
excitons,"*^ and empirical energy-dependent self-energy 
correction according to the Kohn-Sham local density 
theory of excitationji^ However, the simplest method 
is to apply the scissor operator^ which displaces the 
LDA eigenvalues for the unoccupied states by a rigid 
energy shift. Using the latter method excellent agree- 
ment with experiments has been demonstrated for lead 
monochalcogenides^^ and ferroelectric NaN02.^^ How- 
ever, the question as to whether the rigid energy shift 
is generally applicable to semiconductors with strong 
Coulomb correlation effects is open. 

In this work electronic structure and optical properties 
of the ZnX-w and -z phases have been studied in the 
energy range from to 20 eV based on first-principles 
band structures calculations derived from DFT within 
the LDA, GGA, and LDA+U. 



II. COMPUTATIONAL DETAILS 

Experimentally determined lattice parameters have 
been used in the present ab initio calculations (Table H]). 
The ideal positional parameter u for ZnX-w is calcu- 
lated on the assumption of equal nearest-neighbor bond 
lengths)^ 

"-3(c) +4 W 

The values of u for the ideal case agree well with the 
experimental values u* (see Table H]) . Self-consistent cal- 
culations were performed using alOxlOxlO mesh ac- 
cording to the Monkhorst-Pack scheme for the ZnA"-z 
phases, and the T-centered grid for the ZnX-w phases. 



A. Calculations by VASP package 

Optical spectra have been studied based on the band- 
structure data obtained from the VASP-PAW package^^ 
which solves the Kohn-Sham eigenvalues in the frame- 
work of the DFT^^ within LDA^ GGA^ and the sim- 
plified rotationally invariant LDA-|-£/] ^^i^^ The exchange 
and correlation energy per electron have been described 
by the Perdew and Zunger parametrization'^'' of the quan- 
tum Monte Carlo results of Ceperley and Alder. The 
interaction between electrons and atomic cores is de- 
scribed by means of non-norm-conserving pseudopoten- 
tials implemented in the VASP package. The pseudopo- 
tentials are generated in accordance with the projector- 
augmented wave (PAW)^»^ method. The use of the 



PAW pseudopotentials addresses the problem of inade- 
quate description of the wave functions in the core region 
(common to other pseudopotential approaches^) and its 
application allows us to construct orthonormalized all- 
electron- like wave functions for Zn 3d and 4s, and s and 
p valence electrons of the X atoms under consideration. 
LDA and GGA pseudopotentials have been used and the 
completely filled semicore Zn-Sd shell has been consid- 
ered as valence states. 

It is well known that DFT calculations within LDA 
and GGA locate the Zn-3d band inappropriately close 
to the topmost VB, hybridizing with the 0-p band, fal- 
sifies the band dispersion, and reduces the band gap. 
Nowadays the problem is known to be solved by using 
the LDA-f 5/C and LDA-h[/.-2i.S2i26^^^ por the DFT 
calculations within LDA-I-J7 explicit values of the pa- 
rameters U and J are required as input. In previous 
papers^iiS^ we have estimated the values of the U and J 
parameters within the constrained DFT theory^^ and in 
a semiempirical way by performing the calculations for 
different values of U and forcing match to the experimen- 
tally established^^ location of the Zn-Bd bands. Based 
on the result o^^'^^ the values of the parameters U and J 
listed in Table |T] are chosen to study optical spectra. 



B. Calculations by MindLab package 

For investigation of the role of the SO coupling in elec- 
tronic structure and optical properties of ZnX, DFT 
calculations have been performed using the MindLab 
package,— which uses the full potential linear muffin-tin 
orbital (FP LMTO) method. For the core charge density, 
the frozen core approximation is used. The calculations 
are based on LDA with the exchange-correlation poten- 
tial parametrized according to Gunnarsson-Lundquist— 
and Vosko-Wilk-Nussair.^^ The base geometry in this 
computational method consists of a muffin-tin part and 
an interstitial part. The basis set is comprised of linear 
mufhn-tin orbitals. Inside the mufhn-tin spheres the ba- 
sis functions, charge density, and potential are expanded 
in symmetry-adapted spherical harmonic functions to- 
gether with a radial function and a Fourier series in the 
interstitial. 



C. Calculation of optical properties 

From the DFT calculations the imaginary part of the 
dielectric function e2(w) has been derived by summing 
transitions from occupied to unoccupied states for ener- 
gies much larger than those of the phonons: 

'^'^''^ = 2^&^^ J rf'k^(k„|p.|kn')(k„'lp,|kn)x 

nn' 

/kn(l - /kn')^(ekn' - Ckn " (2) 



4 



TABLE I: Theoretically and experimentally (in brackets) determined unit-cell dimensions a, c, volumes V, ideal u (calculated 
by Eq. [T]) and experimental u* as well as values of the parameters U and J from Refs. [U and 112 were used in the present 



calculations. 


For w-type structure a — b. For the 


z-type structure a = b — 


c and all atoms 


are in 


fixed positions. 




Phase 


a, (A) 


c (A) 


V (A^) 


u* 


u 


U (oV) 


J (oV) 


ZnO-w" 


3.244(3.250) 


5.027(5.207) 


y=45. 82(47. 62) 


0.383 


0.380 


9 


1 


ZnS-w' 


3.854(3.811) 


6.305(6.234) 


y=81. 11(78. 41) 


0.375 


0.375 


6 


1 


ZnSc-w*^ 


4.043(3.996) 


6.703(6.626) 


y=94. 88(91. 63) 


0.375 


0.371 


8 


1 


ZnTc-w'' 


4.366(4.320) 


7.176(7.100) 


y=118. 47(114. 75) 


0.375 


0.373 


7 


1 


ZnO-z^ 


4.633(4.620) 




V — yy .-4:0 i^yo.Di J 






o 
o 


■j^ 


ZnS-z^ 


5.451(5.409) 




y=161. 99(158. 25) 






9 


1 


ZnSc-z" 


5.743(5.662) 




y=189. 45(181. 51) 






8 


1 


ZnTc-z'' 


6.187(6.101) 




y=236. 79(227. 09) 






8 


1 
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Here {px,Py,Pz) = p is the momentum operator, /kn 
the Fermi distribution, and |kn) the crystal wave func- 
tion corresponding to the energy ekn with momentum k. 
Since the ZnX-w phases are optically anisotropic, com- 
ponents of the dielectric function corresponding to the 
electric field parallel {E \\ c) and perpendicular {E _L c) 
to the crystallographic c axis have been considered. The 
ZnX-z phases are isotropic, consequently, only one com- 
ponent of the dielectric function has to be analyzed. 

The real part of the dielectric function ei{oj) is cal- 
culated using the Kramer-Kronig transformation. The 
knowledge of both the real and imaginary parts of the 
dielectric tensor allows one to calculate other important 
optical spectra. In this paper we present and analyze 
the reflectivity R{uj), the absorption coefficient a{uj), the 
refractive index n{Lu), and the extinction coefficient k(Lu): 



n{uj) 



k{uj) 



^ei{u) + ei{uj)+ei{uj) 



(3) 



(4) 



(5) 



(6) 



Here e{uj) — ei(w) -1-162(0;) is the complex dielectric func- 
tion. The calculated optical spectra yield unbroadened 
functions, and consequently have more structure than the 
experimental onesi^i^i^i^ To facilitate a comparison 
with the experimental findings, the calculated imaginary 
part of the dielectric function has been broadened. The 
exact form of the broadening function is unknown. How- 
ever, analysis of the available experimentally measured 
optical spectra of ZnX shows that the broadening usu- 



ally increases with increasing excitation energy. Also, 
the instrumental resolution smears out many fine fea- 
tures. These features have been modelled using the life- 
time broadening technique by convoluting the imaginary 
part of the dielectric function with a Lorentzian with a 
full width at half maximum of 0.002(fio;)^ eV and in- 
creasing quadratically with the photon energy. The ex- 
perimental resolution was simulated by broadening the 
final spectra with a Gaussian, where the FWHM is equal 
to 0.08 eV. 



III. RESULTS AND DISCUSSION 

A. Band structure 

The optical spectra are related to band dispersion and 
probabilities of interband optical transitions. So, it is 
of interest to analyze the electronic structure in detail. 
Band dispersions for ZnAT-w and ZnA"-z calculated by 
DFT within LDA and LDA-|-t7 are presented in Fig. [1 
The general features of the band dispersions are in agre- 
ment with previous studies (see, e.g., Refs. [26ll6^ . andpT^). 
It is seen from Fig. [5] that the CB minimum for ZnX- 
w and ZnX-z are much more dispersive than the VB 
maximum, which shows that the holes are much heav- 
ier than the CB electrons in agreement with experimen- 
tal data^^'^^ for the effective masses and calculated with 
FP LMTO,26 LCAO,is as well as our findings.^i'^i Con- 
sequently, mobility of electrons is higher than that of 
holes. Furthermore, these features indicate that p elec- 
trons of X (that form the topmost VB states) are tightly 
bound to their atoms, and make the VB holes less mobile. 
Hence, the contribution of the holes to the conductivity 
is expected to be smaller than that of CB electrons even 
though the concentration of the latter is smaller than that 
of the former. These features emphasize the predominant 
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ionic nature of the chemical bonding. Another interesting 
feature of the band structures is that the VB maximum 
becomes more dispersive with increasing atomic number 
of X from O to Te. 

As noted in our previous contributions, ^^'^^ the band 
gaps of ZnX calculated by DFT within LDA, GGA, 
and LDA+J7 are underestimated and the question as to 
whether it is possible to shift the CB states rigidly was 
kept open. As found from the optical spectra discussed 
on the following sections, rigid shifts of the CB states up 
to the experimentally determined locations can provide a 
good first approximation for the stipulation of the band 
gap. So, for the band dispersions in Fig. [21 we have made 
use of this simple way for correcting the band gaps cal- 
culated by DFT. The only problem in this respect was 
the lack of an experimental band-gap value for ZnO-z. 
To solve this problem, the experimentally and calculated 
(by DFT within LDA) band gaps {Eg) of the ZnX series 
were plotted as a function of the atomic number of X. As 
seen from Fig. [3l Eg for the ZnX-w phases are very close 
to the corresponding values for the ZnX-z phases and 
the shape of the experimental and calculated functional 
dependencies are in conformity. On this basis the band 
gap of ZnO-z is estimated by extrapolating the findings 
for ZnX-z from ZnS-z to ZnO-z. This procedure gave 
Eg « 3.3 eV for ZnO-z. 

It is well-known that not only band gaps are under- 
estimated within LDA and GGA, but also band disper- 
sions come out incorrectly, whereas location of energy 
levels of the Zn-3d electrons are overestimated (see, e.g., 
Refs.lUlllil andlsi). As also seen from Fig. El calcula- 
tions within the LDA-I-C/ approach somewhat correct the 
location of the energy levels of the Zn-3d electrons. The 
elucidation of the eigenvalue problem and the order of 
states at the topmost VB from LDA, GGA, and LDA+JJ 
calculations are discussed in Refs. 20.21,22 and 26 will 
not be repeated here. 

Examination of Fig. [2] shows that the VB comprises 
three regions of bands: first a lower region consists of s 
bands of Zn and X, a higher lying region where well 
localized Zn-3d bands, and on top of this a broader 
band dispersion originating from X-p states hybridized 
with Zn-3c? states. The latter sub-band is more pro- 
nounced in ZnO than in the other ZnX phases consid- 
ered. The hybridization is most severe according to the 
LDA and GGA calculations, whereas the LDA-I-J7 cal- 
culations somehow suppress this and improve the band 
gap underestimation. A more detailed discussion of these 
aspects is found in Refs. [Ill and [2^ 

The SO splitting at the topmost VB is known to play 
an important role for the electronic structure and chem- 
ical bonding of solids. ^^'^^-^^ In semiconductors with z- 
type structure the SO splitting energy is determined as 
the difference between energies of the topmost VB states 
with symmetry Fgi, and Fyi,.— i^Sii^^ In the w-type com- 
pounds the topmost VB is split not only by SO inter- 
action, but also by CF giving rise to three states at the 
Brillouin-zone center. To calculate the SO splitting en- 



ergy for w-type phases quasi-cubic model of Hopfield2^ is 
commonly used. 

It is well known that the SO splitting energy derived 
from ah initio calculations agree well with experimental 
data only for some of the semiconductors. This is demon- 
strated, for example, for all diamond-like group IV and z- 
type group III-V, II- VI, and I- VII semiconductors, w- 
type AIN, GaN, and InN,^ ZnX-w and -z (X=S, Se, and 
Te)^i^ and CdTci^ However, the errors in estimated 
SO and CF splitting energies by LDA calculations are sig- 
nificant for semiconductors with strong Coulomb corre- 
lation effects, as demonstrated, e.g., for ZnOj2ii22ii2& For 
such systems DFT calculations within LDA-t-C/i2LS^i2& are 
shown to provide quite accurate values for Aqf and Aso- 
Overestimation of the p-d hybridization in various vari- 
ants of the DFT can also lead to the wrong spin-orbit 
coupling of the valence bands i^i^l 

Systematic study of the SO coupling parameters was 
performed for zinc-blende II- VI semiconductors (Ref. [30h 
using the TB and LMTO methods, as well as for all di- 
amondlike and zinc-blende semiconductors (Ref. [28l ) us- 
ing the FLAPW method with and without the pi/2 local 
orbitals, and the frozen-core PAW method implemented 
into VASP. The corrections coming out from inclusion of 
the local orbitals are found to be neghgible for the 
compounds with light atoms. Analysis of these results 
shows that the So splitting energy coming out from calcu- 
lations using the VASP-PAW shows good agreement with 
the experimental data. This result was also obtained^i 
recently for ZnX of wurtzite and zinc-blende structures. 
As demonstrated in Refs. [2llj23 the SO splitting energy 
(Ago) increases when one moves from ZnO-z to ZnTe-z 
in agreement with earlier findings of Ref. [28l 

To study the role of the SO couphng in band dispersion 
the present ah initio calculations have been performed 
by VASP and MindLab packages and spin orbit splitting 
energy is found. The results are presented in Table [III 
Analysis of the Table [III shows that (Ago) calculated by 
MindLab is quite accurate. 



TABLE 11: Calculated SO splitting energy (in meV), using 
the MindLab package along with the previous theoretical and 
experimental findings. 



ZnO-z 



ZnS-z 



ZnSc-z 



ZnTc-z 



31 
31 
-34" 
-34' 

-33'' 



66 

66 

66" 

66' 

64" 

64'' 

65" 



432 

432 

393" 

398' 

392° 

393'' 

420-'' 



914 

914 

889° 

916'' 

898"^ 

897'' 

91(/ 



"LAPW, Ref. [23. 
''LAPW-l-pi/2, Ref. [28 
'^VASP-PAW, Ref. ^ 
''VASP-PAW, Ref. "2^. 
'^Experiment, Ref. 73. 
■^Experiment, Ref. 7^. 



As expected, band dispersions calculated with and 
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degenerate states called A, B, and C states with energies 
Eg{A), Eg{B) and Eg{C), respectively)^ located in order 
of decreasing energy, i.e. Eg{A) > Eg{B) > Eg{C). The 
center of gravity of the A, B, and C states, located at 
{Eg{A) — Eg{C))/3 below the topmost A state, remains 
to be nearly the same as the topmost VB corresponding 
to the case without the SO coupling. ^^'^^ Consequently, 
to compare band structures calculated with and with- 
out the SO couphng, one should plot the band structure 
with Fermi energy at the center of gravity of the A, B, 
and C states for the former and at the topmost VB for 
the latter. Hence, when the SO coupling is applied, then 
the A and B states as well as the bottommost CB move 
upwards to {Eg{A) — Eg{C))/2> in energy, whereas the 
C state moves downwards to {Eg{A) — Eg{C))2/3 com- 
pared to the center of gravity. Then, positions of the 
lowest VB region calculated with and without the SO 
coupling remains nearly identical. 



FIG. 3: Band gaps for ZnX-w (circles) and ZnX-z (trian- 
gles) pha ses determined experimentally (filled symbols, from 
Refs. I21I and |23 ) and calculated (open symbols) by DFT 
within LDA as a function of the atomic number of the X 
component of ZnX. 



without the SO coupling differ little when the SO split- 
ting energy is small. However, the difference increases, 
when one moves from ZnO to ZnTe. This feature is 
demonstrated in Table |TT] and Fig. 2] for band disper- 
sions of ZnO-z, ZnO-w, ZnTe-z, and ZnTe-w calculated 
by VASP with and without including the SO coupling. 
As is well known (see e.g. Ref. 21,26,27), without the SO 
coupling top of the VB of ZnX-w is split into a doublet 
and a singlet states. In the band structure, Fermi level is 
located at the topmost one [Fig. [3] , which is the zero en- 
ergy. Upon inclusion of the SO coupling into calculations, 
the doublet and singlet states are split into three twofold 



B. General features of optical spectra of ZnX 

Since optical properties of solids are based on the band 
structure, the nature of the basic peaks in the optical 
spectra can be interpreted in terms of the interband tran- 
sitions responsible for the peaks. Such an interpreta- 
tion is available for semiconductors with z- and w-type 
structures. iiii^i^ In order to simplify the presentation 
of the findings of this work, the labels Eo,Ei, and E2 
of Ref. [ll| (from the reflectivity spectra) were retained 
in Table IIIII and Fig. 21 The subscript is ascribed to 
transitions occurring at F, the subscript 1 to transitions 
at points in the [111] direction, and the subscript 2 to 
transitions at points in the [100] direction (referring to 
the k space for the z-type structure). Assignment of the 
Eq, El, and E2 peaks to optical transitions at high sym- 
metry points is presented in Table IIIII and Fig. [H 
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FIG. 4: Band dispersion for ZnO-z, ZnO-w, ZnTe-z, and ZnTe-w calculated by the VASP-PAW method within LDA accounting 
for SO coupling (solid lines) and without SO coupling (open circles). Topmost VB of the band structure without SO coupling 
and center of gravity of that with SO coupling are set at zero energy. Symmetry labels for some of the high-symmetry points 
are shown for ZnTe-z and ZnTe-w to be used for interpretation of the origin of some the peaks in the optical spectra of ZnX-w 
and ZnX-z. 



The optical spectra ei{uj), e2(w), a{uj), R{uj),n{uj), and 
k{uj) calculated by DFT within LDA, GGA, and LDA+C/ 
are displayed in Figs. [5]-[5] and compared with available 
experimental findings.^'' The spectral profiles are indeed 
very similar to each other. Therefore, we shall only give 
a brief account, mainly focusing on the location of the 
interband optical transitions. The peak structures in 
Figs. [SHH] can be explained from the band structure dis- 
cussed above. 

TABLE III: Relation of the basic Eo,Ei, and E2 peaks in 
the optical spectra of ZnX to high-symmetry points (see 
Refs. I 111 andhj) in the Brillouin zone at which the transi- 
tions seem to occur. 
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All peaks observed by experiments (see, e.g., Refs. [ill 
and Il4l ) are reproduced by the theoretical calculations. 
Because of the underestimation of the optical band gaps 
in the DFT calculations the locations of all the peaks in 
the spectral profiles are consistently shifted toward lower 
energies as compared with the experimentally determined 
spectra. Rigid shift (by the scissor operator) of the op- 
tical spectra has been applied, which somewhat removed 
the discrepancy between the theoretical and experimental 
results. In general, the calculated optical spectra quali- 
tatively agree with the experimental data. In our theo- 
retical calculations the intensity of the major peaks are 
underestimated, while the intensity of some of the shoul- 
ders are overestimated. This result is in good agreement 
with previous theoretical findings (see, e.g., Ref.figf). The 
discrepancies are probably originating from the neglect 
of the Coulomb interaction between free electrons and 
holes (excitons), overestimation of the optical matrix el- 
ements, and local-field and finite-lifetime effects. Fur- 
thermore, for calculations of the imaginary part of the 



dielectric response function, only the optical transitions 
from occupied to unoccupied states with fixed k vector 
are considered. Moreover, the experimental resolution 
smears out many fine features and, as demonstrated in 
Fig. [U there is inconsistency between the experimental 
data measured by the same method and temperature. 
However, as noted in Introduction, by accounting for 
the excitons and Coulomb correlation effects in ab initio 
calculations^ by the linearized-augmented plane-wave 
method plus local orbitals (LAPW+LO) within LDA+U 
not only corrected energy position of the Zn-3d electrons 
and eigenvalues, but also optical response. Consequently, 
accounts for the excitons play an important role in the 
optical spectra. 

The optical spectra calculated within LDA, GGA, and 
LDA+C/ do not differ significantly from each other for 
the ZnX-w and -z phases except for ZnO-w and -z, for 
which the optical spectra calculated within LDA+U are 
significantly different from those obtained by LDA and 
GGA. The difference between the optical spectra calcu- 
lated by LDA and GGA and those calculated by LDA+C/ 
decreases, when one moves from ZnO to ZnTe. For ZnTe 
the difference can be said to be very small. This feature 
shows that in ZnO-w and -z Coulomb correlation effects 
are strong compared to the other ZnX-w and -z phases 
in agreement with recent LAPW+LO calculations^- in- 
cluding electron-hole correlations. 

Comparison of the optical spectra for E ± c and E \\ c 
for each of the ZnX-w phases with the isotropic spectra of 
the corresponding ZnX-z phases shows that the locations 
of the peaks almost coincide. This similarity refiects that 
there is only small differences in the local arrangement 
of the atoms in the ZnX-w and corresponding -z phases. 

C. ZnO-w and ZnO-z 

The optical spectra of ZnO-w and -z calculated by DFT 
within LDA, GGA and LDA+C/, together with measured 
data are displayed in Fig. [5] One clearly sees three major 
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FIG. 5; Optical spectra of ZnO-w for -E || c (first column) and E L c (second column), and ZnO-z (third column). In the 
first and second columns the results obtained from calculations are plotted by thick solid lines for LDA, thinner solid lines for 
GGA, and dashed lines for LDA+(7 as compared to experimental data from Ref. [3 (open circles). In the third column results 
calculated within LDA (solid lines), GGA (dotted lines), and LDA+C/ (dashed lines. a{Lo) is given in cm~^ divided by 10^. ) 




FIG. 6: Optical spectra of ZnS-w for -B || c (first column) and _E ± c (second column), and ZnS-z (third column) calculated 
within LDA (thick solid lines), GGA (thin solid lines), and LDA+(7 (dashed lines) and compared with experimental data (open 
circles) from Ref . [13. a{u)) is given in cm~^ divided by 10^. 
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peaks in the experimental spectra located in the energy 
ranges 3.1-3.3 (Eq), 7.5-8.5 (Ei) and 10-15 eV {E2). In 
the E2 peak £2(1^), ct{uj), R{lu), and k{uj) are seen to take 
larger values than those in the Eq and Ei peaks. This 
is one of the major features, which distinguishes ZnO-w 
and -z from the other ZnX phases. 

It should be noted that with increasing value for the 
parameter U in the LDA+J7 calculations, the intensity of 
the Eq and Ei peaks of ZnO-w decreases compared with 
the LDA and GGA findings as well as with the experi- 
mental data. However, the intensity of the Ei peak of 
ZnO-z from the LDA+C/ calculations has increased and 
has become even larger than those derived from the LDA 
and GGA calculations as well as the experimental data. 
The intensity of E2 peak from the LDA+C/ calculations 
oscillates significantly, showing disagreement with the 
LDA and GGA calculations as well as the experimental 
measurements. Hence, although LDA+U calculations2i 
were good to increase the LDA-derived band gap and the 
SO splitting energy as well as to decrease the crystal- 
field splitting energy and improving the band dispersion, 
it was poorer than LDA and GGA to describe the op- 
tical properties of the ZnO phases. Probably, this dis- 
crepancy comes about because in our ab initio calcu- 
lations electron-hole interactions and SO coupling are 
not includedi^ The strong variation of the optical prop- 
erties with increasing the values of U indicates appre- 
ciable Coulomb correlation effects in ZnO-w and -z in 
agreement with our previous band structure findings2ii22 
and LAPW+LO calculations^*^ including cxcitonic effect. 
This feature is not present in the spectra for the other 
ZnAT-w and -z phases considered. 

For convenience of analysis the e2(cj) profile was ana- 
lyzed by adjusting the peak location to the experimental 
data of Ref. [ij by rigid shift. On comparing this result 
with that of Ref. [sl, it is concluded that the peaks at 
3.40 eV for II c and that at 3.33 eV for E ± c of e2{uj) 
and R{ll!) can be ascribed to transitions at the funda- 
mental absorption edge. As shown in Ref. 8, the energy 
difference (0.07 eV) between these two peaks gives the 
separation between the so-called A, B (for E ± c), and 
C (for E \\ c) states forming the topmost VB of w-type 
semiconductors in agreement with 0.083 eV according to 
the band-structure analyses in Refs. [21I and [2^ 

There are two broad shoulders of the peak Eq located 
at 4.44 and 5.90 eV for E ± c, and 3.90 and 5.29 eV 
for -B II c. Similar shoulders are found at lower energies 
in the experimental spectra of Refs. d and [s^ observed 
at 3.35 and 3.41 eV for ± c and 3.39 and 3.45 eV 
for i? II c, and the origin of these shoulders have been 
ascribed to exciton-phonon coupling. However, in our 
ab initio studies excitons and lattice vibrations are not 
taken into consideration. 



D. ZnS-w and ZnS-z 

The experimental^^ optical spectra for the ZnS-w and 
-z phases are displayed in Fig. [Altogether with those cal- 
culated according to the LDA, GGA, and LDA-|-?7. It 
is seen that the magnitudes of the experimentallyii ob- 
served shoulders around the E2 peak in the reflectivity 
spectra of ZnS-w are overestimated in the DFT calcu- 
lations. As a result, the intensities of the shoulders are 
almost the same as intensities of the peaks Ei and E2 for 
E ± c and even exceeds them for i? || c. 

The calculated optical spectra for ZnS-w by LDA, 
GGA, and LDA-fC/ turned out to be almost identical 
at energies below 10 eV. However, at higher energies the 
LDA- and GGA-derived peaks differ from those obtained 
by LDA+U. This difference can be associated with Zn- 
3d electrons, which were shifted toward lower energies in 
the LDA-t-C/ calculations. Hence, in ZnS-w and -z the 
Coulomb correlation effects appear to play a significant 
role in optical properties at energies higher than 10 eV. 

Compared to ZnO-w, the calculated optical spectra of 
ZnS-w and -z show larger disagreement with the exper- 
imental data. The discrepancy is quite pronounced in 
the absorption and reflectivity spectra of ZnS-w and -z, 
especially at energies exceeding 7 eV. The magnitude of 
the peaks located at higher energies are overestimated 
significantly compared to the experimental data. The 
overestimation is more severe in ZnS-z than in ZnS-w as 
judged from the intensity of the E2 peak. 

E. ZnSe-w, ZnTe-w, ZnSe-z, and ZnSe-z 

The optical spectra for ZnSe-w, ZnSe-z, ZnTe-w, and 
ZnTe-z calculated by DFT within LDA, GGA, and 
LT)A+U are displayed in Figs. [7] and [5] together with 
corresponding experimental spectra. Since experimental 
optical spectra for ZnSe-w and ZnTc-w are not available, 
rigid shift of the parameters toward higher energies have 
been performed on the basis of the reflectivity spectra for 
ZnSe-z and ZnTe-z (Figs. [7] and [8]). A closer inspection 
of Figs. [7] and [5] shows that the optical spectra calculated 
within LDA, GGA, and LDA+U are almost the same for 
all selenide and telluride phases. The small differences 
noted in the absorption and reflectivity spectra appear 
to originate from the Zn-3d electrons. 

The location and magnitude of the experimentally 
measured Ei peak in the reflectivity spectra of ZnSe-z 
and ZnTe-z have been assigned-- '--^^ to fundamental ab- 
sorption and A3 — Ai transitions at the [0.17,0.17,0.17] 
point of the Brillouin zone. These assignments agree 
well with theoretical calculations. However, the theo- 
retical calculations did not locate the Ei + Ai peak on 
the high-energy side of the Ei peak, which was observed 
experimentally for both ZnSe-z and ZnTe-z. The rea- 
son is certainly that SO coupling was not included in the 
calculations. 

The experimental^- Eg peak in the reflectivity spectra. 
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FIG. 7: Optical spectra of ZnSe-w for -E || c (first column) and iJ _L c (second column), and ZnSe-z (third column) calculated 
within LDA (thick solid lines), GGA (thin solid lines), and LDA+(7 (dashed lines) and compared with experimental data (open 
circles) from Ref . [l4l) . a{Lo) is given in cm~^ divided by 10^. 




Energy (eV) 

FIG. 8: Optical spectra of ZnTe-w for £ || c (first column) and i5 _L c (second column), and ZnTe-z (third column) calculated 
within LDA (thick solid lines), GGA (thin solid lines) and LDA+t/ (dashed lines) and compared with experimental data (open 
circles) from Ref. [13. a(a;) is given in cm~^ divided by 10^. 
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corresponding to transitions at k = [viz. from the high- 
est state of VB (Fis) to the lowest state of CB (Fi)] is 
well reproduced by the theoretical calculations. One also 
sees the Eq + Aq peak in the theoretical spectra, which 
previously^'^ was ascribed to SO splitting. Since SO cou- 
pling was neglected in the theoretical calculations, the 
origin of the Eq + Aq peak is not likely to be related to 
SO coupling. 

Similar to the findings for ZnS-z, theoretically calcu- 
lated optical spectra for the lower energy regions of ZnSe- 
z and ZnTe-z agree with experimental findings. How- 
ever, the intensity of the peaks located at higher energies 
are overestimated in the DPT calculations. Anyway, this 
discrepancy is not as severe as that for ZnS-z. The cal- 
culated refiectivity spectra agree well with experimental 
data in the energy range < 6 eV. At higher energies (6- 
15 eV) LDA, GGA, and LBA+U aU overestimate the 
intensity of the reflectivity. Fairly good agreement with 
the experimental data is achieved in the energy range 
15-20 eV for the real and imaginary parts of the dielec- 
tric function for ZnSe-z and ZnTe-z. For the other opti- 
cal spectra of ZnSe-z and ZnTe-z the agreement between 
theory and experiment is poorer. 



F. Influence of spin-orbit splitting on the optical 
spectra of ZnX 

It is well-known that SO splitting at the top of the VB 
of a semiconductor is very important for optical transi- 
tions and one should expect large difference in the optical 
spectra calculated with and without the SO coupling. In 
this section we shall analyze how the SO coupling influ- 
ences the optical spectra of ZnX. For this analysis, ab 
initio band structure calculations have been performed 
using the MindLab software with and without SO cou- 
pling. Based on the band structure studies, dielectric 
response function 62 (w) has been calculated. The results 
for ZnX-z are presented in Fig. [9] and compared with ex- 
perimental data, where it is seen that e2{i^) calculated 
for ZnO-z without SO coupling is slightly higher than 
that with the SO coupling, the main deviations occur- 
ring at 3.44-6.00 eV and 10.00-12.00 eV. The reason for 
the small distinctions in £2(0;) in this case is the small 
SO splitting energy. 

Our findings show that the SO splitting influence the 
calculated optical spectra and in particular it most pro- 
nounced at energies lower than ~ 12 eV. At higher en- 
ergies, the difference between the optical spectra calcu- 
lated with and without SO coupling is fairly small and 
agree reasonable well with the experimental dataJ^ This 
statement applies to all ZnX phases studied. It should be 
noted that intensities of the peaks calculated with the SO 
coupling are generally lower than those obtained without 
SO coupling and the latter set agrees better with experi- 
mental data than the former. Furthermore, in the exper- 
imental spectra there are low intensity peaks located at 
9.4 eV in ZnS-z, 8.4 eV in ZnSe-z, and 7.0 eV in ZnTe- 
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FIG. 9: Imaginary part of the dielectric response function 
e2{Lu) for ZnO-z, ZnS-z, ZnSe-z, and ZnTe-z calculated by 
DFT using the MindLab package within LDA including SO 
coupling into consideration (dashed lines) and without SO 
coupling (continuous lines) and compared to experimental 
data (open circles) from Ref. Q. 



z. However, these peaks are not seen in the calculated 
spectra with the SO coupling. As noted, this discrep- 
ancy can be related to neglect of many above factors like 
Coulomb interaction between electrons and holes, local- 
field effects, and indirect transitions etc. 



G. Role of the ground state structure in the 
optical spectra of ZnX 

In this section we shall analyze optical spectra of ZnX 
calculated using the experimentally and theoretically de- 
termined lattice parameters. To find the lattice param- 
eters from the abinitio calculations, the structural opti- 
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mization has been performed, which includes the follow- 
ing steps: (i) atoms are relaxed keeping the volume and 
shape of the lattice. After convergence is reached, the re- 
sulting lattice and positional parameters have been used 
as input to optimize atomic positions, shape and volume 
of the unit cell altogether. Then dependence of the total 
energy on volume is studied. Minimum of the depen- 
dence was accepted as the equilibrium state. Lattice and 
positional parameters corresponding to the minimum is 
referred to as the theoretically determined lattice param- 
eters. The thus determined theoretical lattice parameters 
do not deviate much from experimental ones. The pa- 
rameters along with the experimentally determined ones 
were used for subsequent computations of the electronic 
structure and optical spectra. The results are presented 
in Fig. [TU]for ZnX-w for E _L. Analysis shows that the 
optical spectra of ZnO-w for E _L deviate each from other 
at lower energies corresponding to the fundamental ab- 
sorption and at higher energies in the range 8-13 eV. 
The reason of the difference can be related to changes 
of the p-d coupling because of the changes of the Zn-0 
bond lengths coming out from structural optimization. 
Optical spectra of ZnO-w for || c and those of other 
ZnX-w and -z calculated using the theoretical and exper- 
imental lattice parameters do not differ each from other 
significantly. 



IV. CONCLUSION 

The band structures of the ZnX-w and -z phases 
(X=0, S, Se, and Te) are calculated by DFT within 
LDA, GGA, and LBA+U. The topmost VB states are 
found to be more dispersive than the bottommost CB 
states. Spin-orbit coupling is found to play an important 
role for band dispersion, location and width of the Zn-3d 
band, and the lowest s band. By analyzing the depen- 
dence of the band gaps on the atomic number of X for 
ZnX, the band gap of ZnO-z is estimated to be ~ 3.3 eV. 
Using the electronic band structures as references, the op- 
tical spectra of ZnX-w and -z are analyzed in the energy 
range 0-20 eV. The locations of the peaks corresponding 
to transitions at the fundamental absorption edge calcu- 
lated by DFT are shifted to lower energies relative to the 
experimental peaks. This deficiency originates from the 
well-known errors in band gaps calculated according to 
DFT. In order to correct the underestimation of band 
gaps calculated by DFT, the location of the calculated 
peaks of the optical spectra have been rigidly shifted to- 
ward higher energies to match with the experimentally 
determined locations. In the thus obtained spectra the lo- 
cations of the peaks in the lower energy region agree well 
with the experimental data. However, the peaks in the 
higher energy region agree only tolerably well with the 
experimental findings. The overall conclusion is that the 
k-independent scissors operator provides a good first ap- 
proximation for correlation of the underestimated band 
gaps for the ZnX-w and -z phases. Based on this result. 
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FIG. 10: Imaginary part of the dielectric response function 
e2{ijj) for ZnX-z calculated by VASP-PAW package using the 
theoretically (solid lines) and experimentally (open circles) 
determined lattice parameters. 



"corrected" band structures of the ZnX phases are ar- 
ranged by adjusting the band gap up to experimentally 
measured value (viz. rigidly lifting the lowest CB). Not 
only the locations, but also the intensities of some of the 
calculated low-energy peaks agree with available experi- 
mental data for all ZnX phases. However, the intensities 
of some peaks located at higher energies and shoulders 
have been overestimated. The GGA approach slightly 
improved the band-gap values. Also, the optical spectra 
of ZnO-w for £' ± c calculated within the GGA agree 
better with the experimental data than those calculated 
within the LDA and LDA-I-J7 approaches. The value for 
the corresponding transition at the fundamental absorp- 
tion edge is decreased and becomes sharper with the use 
of the GGA thus providing better agreement with exper- 
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imcntal data than LDA. For E L c inhomogcncity in the 
electron gas plays an important role, while it is not so im- 
portant for II c. The optical spectra for ZnO-w and -z 
calculated within LDA+JJ for the energy range 0-20 eV 
is found to depend significantly on the location of energy 
levels of the Zn-3rf electrons. For the other ZnX-w and -z 
phases such changes are not so pronounced, in fact, only 
noticeable at energies above 10 eV. Strong Coulomb cor- 
relation effects arc established for ZnO-w and -z. Accord- 
ing to the present LDA-I-J7 calculations the probability 
for the optical transitions at the fundamental absorption 
edge of ZnO-w and -z decrease with increasing U . Optical 
spectra for ZnO-z, ZnSe-w, and ZnTe-w have been pre- 
dicted. The influence of the spin-orbit coupling is found 



to increase with increasing the atomic number of the X 
component of ZnX. 
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